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Abstract 

The Svetitsky-Yaffe conjecture for finite temperature 4cZ SU (2) lattice 
gauge theory is confirmed by observing matching of block spin effective 
actions of the gauge model with those of the 3d Ising model. The 
effective action for the gauge model is defined by blocking the signs of 
the Polyakov loops with the majority rule. To compute it numerically, 
we apply a variant of the IMCRG method of Gupta and Cordery. 
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1 Introduction 



Block spin renormalization group has become an important tool in the 
qualitative and quantitative understanding of critical phenomena in classical 
statistical mechanics and Euclidean quantum field theory. As a basic ingredi- 
ent, it introduces effective Hamiltonians (actions in field theoretic language) 
which govern block spin degrees of freedom. The block spins are determined 
from the original degrees of freedom by averaging them over blocks. 

In principle, renormalization group (RG) solves the problems posed by 
critical or nearly critical statistical systems. Under iterated application of 
the block transformation, either the correlation length in the system becomes 
small, or the flow of effective Hamiltonians eventually reaches a fixed point 
which determines the universal properties of the system. 

A major drawback of the RG approach stems from the fact that effective 
Hamiltonians in general contain an infinite number of couplings, in contrast 
to the original Hamiltonians one starts from, which usually contain only 
a small number of interaction terms. The proliferation of couplings has a 
number of consequences. It is, e.g., not always clear how a certain truncation 
to a finite number of couplings affects the physical results. Furthermore, even 
if one relies on a certain truncation scheme, it might be tedious to explicitly 
compute the effective couplings. 

This might be an explanation why the task of explicit computation of 
block spin effective actions has not received very much attention in the lit- 
erature. See, however, e.g., Refs. p| and |§. 

Svetitsky and Yaffe ||] have conjectured that a (continuous) deconfine- 
ment transition of a (d + l)-dimensional finite temperature lattice gauge 
theory should be in the same universality class as the phase transition of a 
corresponding <i-dimensional spin system. This spin system should have the 
center of the gauge group as a global symmetry group. 

The Svetitsky- Yaffe conjecture for SU(2) gauge theory at finite temper- 
ature offers the possibility of an interesting application of the block spin 
renormalization group. First, it has never been rigorously proved that this 
model belongs to the Ising universality class. On the other hand, the con- 
jecture has been checked several times by comparison of Monte Carlo (MC) 
estimates for the critical indices (which were found in good agreement ||), 
as well as with a mean field like analytical approach (which gives also pre- 
dictions for SU(N) deconfinement temperatures |J). However, so far there 
have been no numerical attempts to explicitly compute the effective action 
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for the Polyakov loops and compare it with that of the Ising model. 

With this article, we intend to fill this gap. We will demonstrate that 
actions for the degrees of freedom relevant for the deconfinement transition 
can well be computed by Monte Carlo. Comparing them with the corre- 
sponding actions for the Ising model we are able to confirm the validity of 
the Svetitsky-Yaffe conjecture in a very fundamental way. 

The article is organized as follows: In Section 2 we introduce the notations 
for finite temperature lattice gauge theory and recall the Svetitsky-Yaffe 
conjecture. In Section 3 we introduce the block spin renormalization group. 
Section 4 explains the idea of flow matching, and the application of Improved 
Monte Carlo Renormalization Group (IMCRG) to SU(2) lattice gauge 
theory is described in some detail. In Section 5 we discuss some details of 
our Monte Carlo methods and present the results. Conclusions follow. 



2 Finite Temperature Lattice Gauge Theory 

Let us briefly review the formulation of finite temperature gauge theory on 
a lattice (see for instance Ref. ||). 

Consider an SU(N) gauge system on a (d + l)-dimensional hypercubic 
lattice of size L d -NT, where L and Nt are the spatial and temporal extensions, 
respectively, in units of the lattice spacing a. 

A Euclidean quantum field theory at finite temperature is obtained if one 
compactifies the (imaginary) temporal direction, keeping infinite the spatial 
directions. In a finite lattice formulation one therefore assumes L ^> Nt- The 
compactification length is proportional to the inverse physical temperature 
T 

N T -a = ^ . (1) 

Denote with U^n) the SU(N) group element belonging to the link with 
origin in the site n = (x, t) and pointing in the /i-direction. The usual 
Wilson action reads 

S g (U) = (3j2(N-neTrU P ) , (2) 
p 

P = ^a<"\ (3) 
9 

where Up is the product of the group elements around the plaquette P. The 
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partition function is given by 



Z 



JldUMexpi-SgiU)} 



(4) 



Because of the periodicity in the temporal direction, the system is also in- 
variant under a global Z N symmetry, i.e. the center of the gauge group: its 
spontaneous symmetry breaking at a finite temperature T c is the signal of 
the deconfinement transition. 

The Polyakov loop is an order parameter for the finite temperature de- 
confinement transition. It is the trace of the ordered product of all timelike 
links with the same space coordinate, wrapping in the time direction 



It is a non-trivial observable from a topological point of view: its vacuum 
expectation value is not invariant under Z^ transformations. It is zero in 
the confining phase, while it acquires a finite value in the deconfined phase. 

According to the 15 years old Svetitsky-Yaffe conjecture |§, integrating 
out the space-like degrees of freedom one should obtain an effective action 
for the Polyakov loops which is short ranged and has the center of SU(N) 
as a global symmetry group. 

Thus, given a (i-dimensional classical spin system with the same sym- 
metry properties, undergoing a continuous phase transition, the (d + 1)- 
dimensional quantum gauge model is expected to be in its universality class 
if the deconfinement transition is a continuous one and the effective Hamil- 
tonian has good locality properties. 

This applies in particular to the 4-dimensional SU(2) gauge model which 
should belong to the 3d Ising universality class. 



To define the block spin transformation, consider a magnetic system consist- 
ing of spins a on the sites of a (i-dimensional lattice, defined by a Hamiltonian 
H and a set of couplings {K}, 




(5) 



t=i 



3 Block Spin Renormalization Group 




(6) 



a 
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The partition function reads 



Z = £exp . ( 7 ) 

M 

The "operators" S a (a) are in general all possible products of spins compatible 
with the symmetry of the Hamiltonian. Explicit examples will be given 
below. 

A block spin transformation maps the fine L d lattice onto the block lattice 
of size L' d , where L = Lb L'. This is achieved by averaging the original spins 
over cubical blocks of side length L B according to a certain rule. 

The Hamiltonian H' of the block spins {/i} assigned to the sites of the 
block lattice is defined by 

exp [—H'(n)\ = J2 P(», <t) exp [-H(a)] , (8) 
M 

where P encodes the mapping from the fine to the coarse lattice. It obeys 
P(/i,ff)>0 and 53P(/x,<t) = 1 . (9) 

The latter property ensures that the partition function remains unchanged, 

Z = £ exp [-#'(//)] . (10) 

M 

In this work we use the majority rule prescription (i.e. the \i spins take values 
plus or minus one) 



(blocks) 

pm= n 2 



1 + /i x / sign °y- 



xex' 



The sign function sign(x) in Eq. (|TTD is defined such that it vanishes for 
x = 0. This ensures that in case of a zero sum of spins inside a block a 
positive (negative) /i x / is selected with probability one half. 

The block Hamiltonian H' can be expressed in terms of operators S' a {fx), 
defined on the block lattice, 

H>(ri = -J2K' a S'M. (12) 
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In the case of the 4-dimensional SU{2) gauge model, the 3d effective ac- 
tion for the signs of the Polyakov loops shares (by definition) the Z 2 symmetry 
with the 3d Ising model. To define this action we assign to each Polyakov 
loop its sign 

a x (U) = sign £(x) . (13) 

Then we block the cr-spins with the majority rule to obtain Ising type block 
/x-spins. 

It follows that, similarly to Eq. (Q), the effective Hamiltonian Ti! for the 
finite temperature gauge system is given by 

exp[-WV)] = J DU P{fji, U) ex P [-S g (U)} , (14) 

with 



(blocks) 



p(^u)= n 2 



(15) 



1 + /v sign (?x(U) 
xex' 

Other procedures of blocking, like first averaging the Polyakov loops inside 
the blocks and then take as Ising type spin its sign, could also be employed. 

We close this section by defining a renormalization group flow. A natural 
way to do it would be to fix a block length, e.g., Lb = 2, and then iterate 
the transformation (|15|). We do not stick to this definition here. Instead 
we define the flow by just increasing the block size Lb- This allows us to 
compute Hamiltonians not only for scales 2 n , but also for arbitrary scales 
Lb, with L B integer. 

4 Monte Carlo Renormalization Group for 
Polyakov Loops 

In this section we first recall the RG matching idea. Then we show how to 
apply the Improved Monte Carlo Renormalization Group (IMCRG) method 
by Gupta and Cordery M to 4d SU(2) gauge theory at finite temperature. 



4.1 Matching of RG Trajectories 

In the infinite-dimensional space of couplings {K}, a renormalization group 
transformation R can be looked at as a mapping of the original bare Hamil- 
tonian H onto a new Hamiltonian H' = R(H), defined by the couplings 
{K'}. 
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Figure 1: Matching of the RG trajectories of two critical models A 
and B in the neighbourhood of an RG fixed point. 



The RG flow obtained under iterated RG transformations will eventually 
end in a fixed point {K*}, defined through H* = R(H*). The critical surface 
is identified by all flows connected to the fixed point in this way. 

The RG matching method is based on the assumption that different phys- 
ical systems, belonging to the same universality class, will follow RG flows 
which originate from different "bare couplings" on the critical surface and 
eventually match in a neighbourhood of the common fixed point. Of course, 
a matching close to a non-trivial fixed point will only take place if both 
models under consideration are at criticality. A matching thus confirms both 
universality and allows to check for criticality. This matching method has 
been successfully applied in the context of spin models, see e.g. 0. The fea- 
sibility of matching different critical flows by means of MC methods mainly 
relies on the assumption that the different trajectories come close to each 
other, i.e. approximately match, before the fixed point is actually reached, 
c.f. Figure [TJ. As we shall see, this condition is met for our models. 

4.2 Effective Couplings from IMCRG 

The Improved Monte Carlo Renormalization Group method ]7| allows to 
compute effective actions for Ising type block spins.[] 

1 A generalization to systems with continuous block variables is not straightforward. 



6 



The main idea is to avoid simulations of the original partition function. 
Instead, consider a modified system defined through 



where 



M 

= EE P ^°) exp[-H(a) + H(n) 



(16) 



(17) 



is a guess for H'(/i). 

This system can be simulated once H and H are given. Note the plus 
sign in front of H in Eq. (|16|). 

The system with partition function Z c is non-critical, even in the case of 
a critical Hamiltonian H(a). Consider the expectation values 



<S' a>c 



-fT'(/i) + 



1 



-H(a) + #(/i) 



(18) 
(19) 



If the guess is exact, i.e., 

#0*) = #'(//) , (20) 

the block spins /i x / completely decouple and fluctuate independently. In other 
words, the system is non-critical and the correlations in Z c are bounded by 
the block size Lb- The correlations functions are then known exactly, 



< S a > 
< S a Sp > 



o, 

n a Sap , 



(21) 



where trivial multiplicity factors. 

Let us assume that H(n) is close to H'(fx). Then a first order expansion 
gives 



< S' a > c = n a (K' a - K' a ) + 0((K' a -K t 



W \2 



(22) 

Solving this equation for K' a allows to improve the guess K' a . Usually a few 
iterations 



K' + n- 1 < S' > c 



(23) 
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where the expectation values are determined by simulation of the system ( |TED , 
are sufficient to determine H' to a good precision. 

To apply the IMCRG procedure to the SU(2) gauge system, one has to 
simulate the partition function 



£ exp [-H'(fi) + HQi.) 
M 

Y,I DU P ^ U ) ex P [S g {U) + H(fl) 

M 



(24) 



Remember that the /i-variables are defined as the majority rule block spins 
of the signs of the Polyakov loops. It is straightforward to design a MC 
procedure for the updating of this system. Details will be given in the next 
section. 



5 Monte Carlo Simulations 

We applied the IMCRG method to three different systems: the 3d standard 
Ising model (with nearest-neighbour coupling), the 3d "J3" model, which 
includes also a third, cube-diagonal neighbour coupling fl(|, and the Ad 



SU(2) pure gauge model. We simulated the system defined by Eq. QIE| ) for 
the spin models and by Eq. ([24]) for the SU (2) gauge model. 

For the updating of the Ising model we used a Metropolis algorithm: A 
single spin a x is proposed to be flipped. It is checked whether this update 
leads to a flip of the block spin /i x /, with x e x'. The total change of 
energy AH (a) — AH(fj,) is then computed and used in the usual Metropolis 
acceptance/rejectance step. 

In case of the SU (2) model, only the temporal links couple to the com- 
pensating block Hamiltonian. The space-like links are updated using the 
incomplete Kennedy-Pendleton heat bath sweep [Tl[ supplied with a number 



of overrelaxation sweeps. For temporal links one employs again a Metropo- 
lis procedure: A proposed change of a link matrix leads to a change of the 
Polyakov loop £(x) of which it is member. If the sign cr x changes, the block 
spin /i x / in turn might flip and give rise to a change of H(fi). The relevant 
energy change for the Metropolis step is AS g (U) — AH{y). 

In practical calculations one has to truncate the interactions in H' and 
H . We chose to include in the ansatz eight 2-point couplings and six 4-point 
couplings. The 2-point couplings can be labelled by specifying the relative 
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Kl2 Kl3 Kl4 

Figure 2: Graphical definition of 4-spin couplings included in the 
effective Hamiltonian. 



position of the interacting spins (up to obvious symmetries): Our couplings 
^...Ks then correspond to 001, Oil, 111, 002, 012, 112, 022, 122. The 
4-point couplings Kg . . . K u are defined in an obvious way through Figure ||. 
The corresponding interaction terms in the effective Hamiltonian are denoted 
byS' a ,a = l...U. 

5.1 Reduction of Critical Slowing Down 

A merit of the IMCRG method is that block spin observables are (nearly) 
decorrelated and the critical slowing down problem is less severe than in 
standard simulations. 

In Figure ^ we show scatterplots (MC time history) of measurements of 
the nearest neighbour block spin 2-point function. The comparison is be- 
tween a simulation of the pure gauge system (without IMCRG compensation 
on the block level) and the system defined through Eq. (^). The plot clearly 
shows that the IMCRG type simulation suffers much less from critical slowing 
down. 

Analogous observations were made in case of the Ising model simulations 
with the compensating block Hamiltonian switched on. It is the reduction 
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Figure 3: Comparing scatter plots of the nearest neighbour block 
correlation function observable in SU(2) simulations with Nt = 2. 
The standard simulations (no compensation on block level) are shown 
on the left, the IMCRG simulations on the right. 
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model 


couplings 


L B 


standard Ising 


0.2216544 


3,5,7,9,13 


I 3 Ising 


(0.128003, 0.051201) 


3,5,7,9,13,17 



Table 1: Block sizes Lb and couplings used for the IMCRG sim- 
ulations of the standard and I3 Ising models. For all runs we used 
L' = 8. 



of critical slowing down obtained from the compensation even on moderate 
lattice sizes which enabled us to obtain reasonable results with moderate 
CPU expense. 

5.2 Matching of the Two Ising Models 

We started by comparing the RG flow of the two Ising models. In Ta- 
ble H we summarize some parameters of the MC simulations. We made 
simulations on lattices consisting of 8 3 blocks of size at the infinite 



volume critical couplings (3 = 0.2216544 ]nj for the ordinary Ising and 
((3 u f3 3 ) = (0.128003,0.051201) for the J 3 model. At each RG step (fixed 
Lb value) we made usually two, three or four IMCRG iterations in order to 
have the guesses of the effective couplings converge to reasonable precision. 
The number of sweeps in each run ranged from ~ 10 5 for the largest lattice 
sizes to a ~ 6 ■ 10 6 for the smallest ones. 

As the final error of the estimate for an effective coupling, we took the 
maximum of the statistical error and the last change of guess in the IMCRG 
procedure. In Figure [| we show the results for the flow of the two leading 
couplings K[, and K' 2 , with increasing block size Lg. To achieve matching, 
the block sizes L B of the J 3 model have been rescaled by a factor of A = 0.59. 
A rescaling is always needed to obtain matching. The reason for this is that 
the two models have a different distance to "travel" before they meet on a 
common trajectory, c.f. Figure |l[ 

The figure shows that the two flows collapse nicely on a single trajectory, 
indicating that they are approaching a common fixed point. This happens 
also with the other 12 couplings not shown in the plot. It turns out that the 
approach to the fixed point can be well fitted by a power law, 

K a (L B ) = K* a + a a -L B p . (25) 
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flow of K_l 




2 4 6 8 10 12 14 
L_B 
flow of K_2 




Figure 4: Flows of nearest and second nearest neighbour couplings in 
the standard (diamonds) and the ^3 (bars) Ising model with increasing 
block size Lb- In order to obtain matching, the block sizes of the I3 
model were rescaled by a factor A = 0.59. The dotted lines are fits of 
the flows with a power law. 
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a, 


model 


L B > 


K* 

a 




P 


X 2 /dof 


001, 


Ising 


3 


0.1990(5) 


-0.077(4) 


1.37(7) 


0.99 


001, 


Ising 


5 


0.1977(7) 


-0.15(7) 


1.9(3) 


0.01 


001, 


h 


3 


0.1981(4) 


-0.240(8) 


1.60(4) 


0.50 


001, 


h 


5 


0.1975(5) 


-0.33(9) 


1.80(14) 


0.09 


001, 


combined 


1 


0.1979(b) 


—0.10(3) 


l.b7(20) 


0.17 
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-0.27(9) 






Oil, 


Ising 


5 


0.0224(1) 


0.023(2) 


1.67 fix 


0.94 


Oil, 


h 


5 


0.0225(2) 


0.055(6) 


1.67 fix 


0.27 


111, 


Ising 


5 


0.0013(1) 


0.015(2) 


1.67 fix 


0.15 


111, 


h 


5 


0.0013(1) 


0.029(2) 


1.67 fix 


1.12 


002, 


Ising 


5 


-0.0202(1) 


0.047(2) 


1.67 fix 


1.53 


002, 


h 


5 


-0.0201(3) 


0.059(8) 


1.67 fix 


0.38 


9, 


Ising 


5 


0.00210(4) 


-0.006(1) 


1.67 fix 


0.10 


9, 


h 


5 


0.00215(10) 


-0.005(3) 


1.67 fix 


0.17 



Table 2: Fit results for the Ising model flows for a number of 2-point 
couplings and for the largest 4-point coupling Kg. The fits were done 
with Eq. (p5[). The second column gives the minimum block sizes that 
were used in the fit. K* are the estimates for the fixed point values. 
A "fix" after a parameter means that the value was kept fixed during 
the fitting procedure. 



If the flows of two models obey such a law, with the same fixed point K* a 
and exponent p, but different "amplitudes" a and a', the rescale factor A to 
obtain matching is 




A= - . (26) 



We always used as a reference trajectory the flow of the standard Ising model 
and rescaled the block sizes of the other models by an appropriate factor. 

The results of our various power law fits of the Ising model are summa- 
rized in Table |2|. We fitted the models separately, checking also the effect of 
discarding the effective couplings for the smallest block sizes. That the fixed 
point value and the exponent of the two models coincide is confirmed by a 
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N T (3 L B L 3 -N t 



1 0.8730 2,3,4,5,6,7 42 3 

2 1.871 3,4,5,6 36 3 
2 1.874 3,4,5,6 36 3 
2 1.877 2,3,4,5,6 36 3 
2 1.880 2,3,4,5,6 36 3 



1 

2 
2 
2 
2 



Table 3: Lattice sizes and values of (3 used for the 4d SU(2) model 
at finite temperature. In the last column the maximum lattice size is 
shown. For all runs we used L' = 6. 



common fit of the two Ising flows, where in the fit function only the ampli- 
tudes a a were allowed to depend on the model. This yields the results quoted 
in the last two lines of the first block of the table. We then fixed p = 1.67 
and fitted the flows of the non-leading couplings with two parameters (fixed 
point value and amplitude of power law correction). We found a very nice 
agreement of the resulting fixed point values for all couplings. 

The value of the exponent p turns out to be too big to be identified with 
the first correction to scaling exponent u ~ 0.8. We expected that to should 
be the leading exponent. A possible explanation of the present observation is 
the following: The amplitude of a power term with u as exponent is too small 
to be detected within our precision. The exponent p with its relatively large 
amplitudes is due to the presence of a redundant operator of the particular 
blocking scheme we used. 



5.3 Matching of the 4d SU(2) with the Ising Models 



We then turned to the 4c? SU(2) gauge model at finite temperature. Infor- 
mations on MC simulations made in this case are given in Table |||. We made 
MC simulations for the Nt = 1 and Nt = 2 cases on lattices consisting of 6 3 
blocks of size L B . 

As the critical deconfinement transition value we used the gauge couplings 
(3 C = 0.8730(2) m for N T = 1. For N T = 2 we studied a neighbourhood of 
the critical value (3 C = 1.880(3) [0 (see Table |). 

For the SU (2) model statistic has necessarily been reduced compared to 
the Ising models. For Nt = 2 measurements ranged from ~ 10 4 for the 
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0.195 




0.155 



Figure 5: Flow of the nearest neighbour coupling in the effective 
action for Nt = 1 lattice gauge theory with gauge coupling (3 = 0.8730 
(squares). Also shown is the standard Ising (diamonds) and and I3 
model (bars). The rescaling factor of the gauge block size with respect 
to the standard Ising scale is A = 0.61. 



largest sizes up to ~ 5 • 10 5 for the smaller ones. 

The K[ leading coupling result for Nt = 1 is shown in Figure |5|. Also 
shown is the Ising flow. Notice that here and in the figures which follow the 
Ising flows are the same as in Figure |], which also means that the fit lines 
plotted are those obtained from the Ising data. Fits were made again ac- 
cording to Eq. ( p5|) . Omitting the L B = 2 value, they gave results consistent 
with the fits of the Ising data, both for the exponent and the asymptotic 
values, however with bigger errors due the lower statistics. The rescaling of 
the gauge block sizes with respect to the standard Ising scale used in this 
case is A = 0.61. 
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-0.005 



-0.01 - 



-0.02 



-0.025 







+ 




i *t 









2 4 6 8 10 12 14 
L B 




Figure 6: Flows of four different effective couplings for Nt = 2 
SU{2) lattice gauge theory at (3 = 1.877, matching with the two Ising 
models (bars, diamonds, and fit lines). The block sizes of the gauge 
model are rescaled by a factor A = 0.65 with respect to the standard 
Ising scale. 
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a Lb = 3 Lb = 4 Lb = 5 Lb = 6 



001 


15962(941 

KJ • IUl/UiJI iJ ~ y 


[501 


17442(52) 

KJ • J. J A. £u \ KJ ^ I 


[341 


18333(70) 


[271 


18800(92) 

KJ • 1UUUU 1 l/j— J ) 


[481 


Oil 


02918(57) 


[1181 
L J 


02800(32) 

V / • \J U KJ V 7 V J \ KJ I 


[41 


02689(42) 


[251 


02500(57) 

KJ > KJ U KJ KJ KJ \ KJ 1 / 


[581 


111 


0.00668(73) 


[1031 

L J 


0.00432(35) 


[24] 


0.00357(45) 


[68] 

L J 


0.00300(61) 


[801 

L 1 


002 


-0.00728(94) 


[177] 


-0.01282(41) 


[62] 


-0.01514(49) 


[77] 


-0.01691(69) 


[139] 


012 


-0.00277(46) 


[70] 


-0.00444(21) 


[12] 


-0.00500(25) 


[56] 


-0.00505(38) 


[35] 


112 


-0.00115(42) 


[62] 


-0.00177(21) 


[5] 


-0.00179(26) 


[59] 


-0.00164(35) 


[53] 


022 


0.00042(63) 


[131] 


-0.00041(29) 


[59] 


-0.00041(39) 


[14] 


-0.00077(49) 


[77] 


122 


-0.00030(52) 


[33] 


0.00012(21) 


[42] 


-0.00043(30) 


[55] 


0.00046(36) 


[110] 


9 


0.00085(34) 


[31] 


0.00146(17) 


[20] 


0.00182(21) 


[48] 


0.00212(28) 


[45] 


10 


0.00005(20) 


[13] 


0.00014(10) 


[20] 


0.00014(12) 


[11] 


0.00020(16) 


[3] 


11 


0.00026(55) 


[32] 


-0.00096(28) 


[70] 


0.00006(33) 


[23] 


-0.00040(43) 


[8] 


12 


-0.00013(18) 


[6] 


-0.00013(10) 


[9] 


-0.00013(12) 


[11] 


-0.00030(16) 


[10] 


13 


-0.00020(36) 


[9] 


-0.00020(20) 


[15] 


-0.00029(23) 


[41] 


0.00004(31) 


[66] 


14 


0.00012(59) 


[4] 


0.00012(33) 


[5] 


0.00012(39) 


[11] 


0.00033(53) 


[95] 


stat 


12 • 10 3 




27 • 10 3 




25 • 10 3 




20 • 10 3 





Table 4: Values of the effective couplings for the 4c? SU(2) model at 
Nt = 2, (3 = 1.877 for different block size Lb- In the bottom row the 
statistics of the last IMCRG iteration (fixed Lb) is given. Statistical 
errors are given in parenthesis. Square brackets contain the change of 
the coupling AK' a in the last IMCRG iteration. 



The flows of the effective couplings for Nt = 2 at j3 = 1.877 are given 
in Figure |6|. One can see a clear matching with the two Ising models (bars 
and diamonds) and the fit lines for the first two couplings. The L#'s of the 
SU(2) gauge model were rescaled in this case by a factor of A = 0.65. 

We found that, among the different gauge couplings used for Nt = 2, 
the supposed critical coupling (3 = 1.880 is actually ruled out, no matter 
which value of the rescaling parameter A is chosen. This can be clearly seen 
in Figure [7] where the flows of the NN coupling in the effective action are 
shown for the four different gauge couplings (3 = 1.880, 1.877, 1.874, 1871 and 
compared with the fitted curve of the Ising model. 

For (3 = 1.880 the system is definitely in the deconfined phase, whereas 
the flow for (3 = 1.871 moves away towards the high temperature fixed point 
(confinement phase). The effective coupling values for the best matching 
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K_l 

0.2 
0.19 
0.18 
0.17 
0.16 
0.15 
0.14 

1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 
L_B 

Figure 7: Flows of the NN coupling in the effective action for Nt = 2 
4d SU{2) at four different gauge couplings (3 = 1.880, 1.877, 1.874, 1871 
(triangles, diagonal crosses, squares and stars respectively). Also 
shown is the Ising flow (diamonds, bars and fitted curve). The rescal- 
ing of the gauge block sizes with respect to the standard Ising scale is 
A = 0.65. 




18 



a 


V = 8 


V = 6 


001 


0.19532(20) 


0.19529(60) 


Oil 


0.02307(6) 


0.02315(13) 


111 


0.00172(14) 


0.00192(51) 


002 


-0.01899(7) 


-0.01946(16) 


012 


-0.00551(5) 


-0.00558(11) 


112 


-0.00168(3) 


-0.00168(11) 


022 


-0.00002(19) 


0.00002(13) 


122 


0.00010(3) 


0.00024(29) 


K' g 


0.00195(2) 


0.00191(9) 


K' w 


0.00024(2) 


0.00024(4) 


K' n 


-0.00050(3) 


-0.00049(11) 


K'l2 


-0.00010(1) 


-0.00007(9) 




-0.00014(8) 


-0.00010(8) 


K'u 


-0.00003(7) 


-0.00025(29) 


stat 


1.3 • 10 6 


3- 10 5 



Table 5: Comparison of effective critical couplings for different sizes 
of the coarse lattice. The example shown is the ordinary Ising, with 
Lb = 9. In the bottom line the statistics is given. 



trajectory of (3 = 1.877 are given in Table f|. The less significant values, 
at Lb = 2, have been omitted. We explicitly reported the statistical errors 
and the AK' a variations in the last IMCRG iteration (the latter in square 
brackets) . 

Let us conclude this analysis with two remarks. 

First, it is worthwhile to stress that weak finite size effects are present 
within this approach: Comparing the 6 3 block lattice of the gauge model with 
the 8 3 of the Ising model should not give sizable systematic errors within our 
precision. 

As a check, in Table || the effective coupling values of the ordinary Ising 
model are reported for two different block lattice sizes, V = 8 and V = 6. 
The result confirms that all couplings are consistent within errors. 

Finally, let us notice that within our statistic the (3 = 1.874 flow can 
also be made compatible with the Ising trajectory: A better resolution to 
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discriminate between the two beta values would have required to extend 
the MC analysis to bigger block sizes Lb, of course with much more CPU 
time consuming. Even though, using the block sizes at our disposal the 
corresponding fit is not as good as that of the (3 = 1.877 value. 

Therefore we assume the latter as the critical coupling value for Nt = 2, 



consistently (within one standard deviation) with Ref. [14 



6 Conclusion and Outlook 

The discussion of MC results shows that the Svetitsky-Yaffe conjecture is 
confirmed in a very fundamental way by observing matching of the SU(2) 
RG trajectory with that of the Ising model. 

At the same time, we showed that IMCRG works well as a method to 
compute the effective action of Ising type degrees of freedom in a genuine 
non-Ising model like Ad finite temperature SU(2) gauge theory. 

Notice also that this kind of calculations could be done on workstations, 
with relatively small computer resources. 

An extension to Nt greater than two would be interesting but more ex- 
pensive. The reason is that with increasing temporal size the small Lb actions 
move farther away from the fixed points, i.e. they need to be blocked more 
in order to come close to the reference Ising flows. This observation is in 
agreement with the fact that also in more standard approaches, e.g. via the 
Binder cumulant, the spatial size of the lattice has to be increased very much 
with increasing Nt- 

Finally, it would be of interest to check this approach with different block- 
ing prescriptions. The rate of approaching the RG fixed point is in fact very 
sensitive to the blocking rule used and a faster convergence can in principle 
be obtained using a more sophisticated blocking scheme than the majority 
rule. 
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